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MAGNETOSPHERIC PHENOMENA 


Juan G. Roederer 


I. INTRODUCTION AND WARNING 

This is supposed to be a Rapporteur Paper, reviewing the contri- 
butions on Physics of the Magnetosphere presented at the 9th IUPAP 
International Conference on Cosmic Rays, and placing them in the 
light of present knowledge in this field. 

This, however, is a very difficult task to accomplish properly. 
First of all, the study of magnetospheric phenomena is only indirectly 
related to what one now has adopted de facto as "Cosmic Ray Phy- 
sics. M This explains the comparatively small number of papers pre- 
sented at this meeting on the subject, and the quite unbalanced n nat- 
ional distribution" of authors ( some countrie s, very active in this field, 
deliberately did not submit contributed papers on this subject). Second, 
the field itself is so wide, that it never could have been covered prop- 
erly in the frame of only a subsection of a conference. Finally, this 
field has progressed so much in the last few years, that a rapporteur 
paper, if given in the traditional way before a cosmic ray audience, 
would require a long introduction just to explain some of the new basic 
concepts on which recent research is centered. 

All these considerations constitute a reasonable excuse for depart- 
ing from tradition in the presentation of this rapporteur paper: (a) 
contributed papers presented at this conference will not be mentioned 
explicitly, neither by title nor by author; (b) no specific references 
will be given at all; (c) only a qualitative discussion will be presented, 
describing the physical processes involved, rather than listing quanti- 
tative experimental and theoretical results; (d) no fancy figures will be 
reproduced; only qualitative sketches will be given, whenever necessary. 
For all these reasons, this paper should never be quoted for reference. 

Two excellent, review papers have been given at this Conference by 
N. F. Ness 1 and S. N. Vernov. 2 In addition, several chapters in two 
recent books, edited by H. Odishaw 3 and W. N. Hess, 4 respectively, 
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deal in a detailed and quite up to date way with this subject, containing 
abundant literature references. For a quite complete collection of ref- 
erences, consult. 5 Many of the following discussions are based on the 
lectures given at the NATO Advanced Study Institute on Radiation Trap- 
ped in the Geomagnetic Field, held in Bergen, Norway (August 1965). 6 

In the qualitative description which follows, we shall start out in 
the solar wind, penetrate the stationary shock wave into the transition 
region, cross the magnetopause, work our way through the distorted 
field in the outer magnetosphere, diffuse into the trapping region and 
to lower altitudes, to finally end up precipitating into the upper atmos- 
phere . 

II. SHAPE AND STRUCTURE OF THE MAGNETOSPHERE 

Let us start out in the solar wind. We now know that there is an 
uninterrupted flow of plasma radially outwards from the sun, repre- 
senting the continuous expansion of the corona caused by a steady dis- 
sipation of energy into heat at its base. Recent satellite measurements 
have established quite convincingly that the solar wind velocity never 
decreases below a certain "base line" of about 320 Km / sec, not even 

during very quiet periods at solar 
minimum. Transient increase s of 
the velocity occur in correlation 
with active regions on the sun; after 
their passage, the plasma flow falls 
down to its base line very readily. 
Ion fluxes are of the order of 3 x 10 8 
cm“ 2 sec' 1 , with kinetic tempera- 
tures of the order of 10 5o K at 1 a.u. 
Ion composition is mainly hydrogen; 
helium ions have been recently 
identified. 


Embedded in this plasma, 
there is the familiar spiralled radial magnetic field, frozen in and 
carried outwards by the infinitely conducting medium. During quiet 
conditions, the field intensity is about 5 gammas. In correlation with 
the passage of active regions on the solar disk, this intensity may 
increase to 100 gammas. The average time lag of 4.5 days between 
the central passage of the active region and the field increase at 1 a.u. 
is in good agreement with the direct measurements of plasma velocity, 
and definitely proves the solar origin of the interplanetary field. 
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Under these conditions of density, temperature, field and bulk 
velocity, the solar wind represents a collisionless plasma flowing at 
supersonic speed with Alfven Mach number 6-7 during quiet periods. 
To this flow, the earth’s magnetic field is interposed as an obstacle. 

The first n warning n given to a solar wind particle, that an obsta- 
cle is being approached, comes in form of a stationary shock wave . In 
ordinary gas flows at supersonic speed, the information about an obsta- 
cle is transmitted to incoming fluid particles by collisions with fellow 
particles in a very thin shock layer, where non-adiabatic compression 
and heating, and a general slowing down to subsonic velocities take 
place. This region is typically a few collision mean free paths thick. In 
a collisionless plasma, information can only be transmitted by the mag- 
netic field, which here replaces the collision process in its key role of 
tying particles together . In a collisionless shock, a similar compression 
and non-adiabatic heating is likely to occur, in a region about one ion 
cyclotron radius thick. This would be, for typical solar wind parameter s, 
of only 1000 Km. 

The existence of this shock or bow wave is now perfectly estab- 
lished experimentally, and its position and shape in the "front* 1 side of 
the magnetosphere (and up to about ± 120° of the earth-sun line) are quite 
well determined near the equatorial plane. The geocentric distance to 
the sub solar point of the shock is about 13 earth radii during quiet con- 
ditions. The experimental study of this shock is of considerable value 
for plasma physics, since it provides essential information for the yet 
unsolved theoretical problem of collisionless shocks. 

Right behind the shock, satel- 
lites and space probes indeed have 
identified a transition region of sub- 
sonic plasma flow, heated up to tem- 
peratures of a few million degrees. 

An important expected effect of the 
shock is that it will compress, i.e. 
amplify in an extremely short time 
any small irregularity of the inter- 
planetary field flowing transversally 
into it; these irregularities may 
then contribute to a general break- 
up into turbulence right behind the 
shock, in the transition region. This 
turbulence is badly needed as the 
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randomizing process which can account for the entropy increase during 
the non-adiabatic compression of the collisionless gas. 

Experimental evidence for turbulence in the transition region was 
found in magnetic field measurements, which show chaotic changes in 
field direction and intensity. On the other hand, plasma measurements 
reveal a sudden change from a directional, steady flow into a consider- 
ably isotropized, irregular flow, whenever the shock is traversed from 
interplanetary space into the transition region. 

The turbulence in the transition region may be responsible for 
statistical (Fermi-type) acceleration of electrons and protons. On the 
other hand, the sudden compression of field irregularities in the shock 
may provide a betatron mechanism for particle acceleration. Evidence 
for sporadic occurrence of intense fluxes of >40 kev electrons was found 
in the transition region. On the other hand, there is evidence for a non- 
Maxwellian high energy tail in the proton distribution. These effects are 
stronger at the dawn side of the transition region. 

The transition region acts like an elastic medium transmitting the 
kinetic pressure of the solar wind right onto the geomagnetic field: the 
infinitely conducting plasma in the transition region will push and com- 
press this field right up to a point where balance is achieved between 
the kinetic pressure and the geomagnetic field pressure. Again, the 
thickness of the boundary layer will be given by a typical ion cyclotron 
radius in the compressed geomagnetic field. Impinging particles from 
the thermalized solar wind in the transition region will tend to drift for 
some time along the boundary, originating electric currents. These are 
precisely the currents necessary to confine the earth f s field in a finite 
volume. The sudden termination of the geomagnetic field is called the 
magnetopause . 

There is now ample experimental information about the magneto- 
pause, its thickness, position and shape near the equatorial plane, in 
the front side of the magnetosphere. Field measurements suggest values 
as low as 100 Km for the boundary thickness; a typical position of the 
sub solar point of the magnetopause is 10 earth radii, during quiet times. 
So far, there is good agreement between these measurements and recent 
theoretical models for the boundary, at least near the equatorial plane. 
These models try to determine the shape and position of a surface such 
that the (known) kinetic pressure of the plasma outside equals to the 
magnetic pressure inside, given by the (known) earth’s internal field, 
plus the field caused by the currents flowing in the surface, which must 
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be placed there in such a way as to cancel exactly the total field outside . 
Once the position of the magnetopause is known, one can compute the 
position of the shockwave, and herewith test the various assumptions 
made in the theoretical treatment of the collisionless shock. 

So far there is no experimental information about the magneto- 
pause at high latitudes, in the noon meridian plane. This, however, is 
a very interesting region, since it is there, where neutral points in the 
magnetic field are to be expected. These neutral points are of consid- 
erable importance for the mechanisms of particle transfer from the so- 
lar wind to the magnetosphere. 

The magnetopause extends into the anti- sun direction encircling 
the tail of the magnetosphere; measurements were made up to about 3 0 
earth radii. The plasma motion in the transition region far away from 
the stagnation point seems to be again highly directional, flowing along 
the magnetopause. There is good evidence for a 5° tilt of the whole tail 
with respect to the earth- sun line, probably related to the interplane- 
tary magnetic field, which, due to its twist, increases the total pres- 
sure on the magnetosphere at the dawn side. 

An extremely important question is that of the stability of the mag- 
netopause. It is now believed that instabilities at the boundary are the 
origin of the geomagnetic activity, i.e. the K p index. They probably 
cause h.m. waves to propagate down to the earth, showing up as fluctu - 
ations in the surface field. Essentially three different situations have to 
be analyzed: quiet conditions (stationary boundary), increase of solar 
wind (compression), and decrease of solar wind (expansion to the normal 
state). During absolutely quiet conditions in interplanetary space, the 
turbulence behind the shock may be 
a persisting cause for instabilities, 
giving rise to the quiet time rema- 
nent geomagnetic activity, especial- 
ly in the polar caps. Stability con- 
ditions during compres sion could be 
entirely different from tho se holding 
during expansion. In effect, recent 
solar wind measurements seem to 
show a correlation of A p index in- 
creases with the increasing phase 
of the solar wind flow, and not with 
the absolute value of the solar wind 
velocity. 
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Instabilities in the magnetopause must play a crucial role in two 
important processes: energetic particle transfer across the boundary, 
and viscous -type momentum transfer from the solar wind to the mag- 
netosphere . 

We now come to the field structure inside the geomagnetic cavity . 
The sources of this field are: the magnetization of the earth 1 s interior, 
the currents flowing on the surface of the magnetopause, the currents 
in the neutral sheet in the tail of the magnetosphere (see below), and, 
eventually, diamagnetic ring currents originating in trapped particle 
density gradients at 2-4 earth radii, in the equatorial plane. During 
quiet conditions, the latter may be of very small effect. At geocentric 
distances less than, say, 5 earth radii, only the internal geomagnetic 
field dominates; beyond 5 R e , the currents in the magnetopause (and in 
the neutral sheet) perturb the dipole -type internal field, and introduce 
a strong noon-midnight asymmetry. 

Most remarkable is the recently established radial character of 
the field in the tail of the magnetosphere, and the existence of a thin 
layer or "neutral sheet", which separates the two regions of mutually 
opposite field directions. The orientation of the neutral sheet is con- 
trolled by both, the solar wind and the earth’s dipole axis: it wobbles 
about the sun-earth line, keeping its normal in the sun-earth-dipole - 
axis plane. 


The lines of force stretching 
out into the tail, come from the po- 
lar regions of the earth. Consequent 
ly, there are two quite distinct fam- 
ilies of field lines emerging out of 
the earth’s surface: those which re- 
turn back to earth, and those which 
are "lost" into the tail. No particles 
can ever be trapped on the latter. 
Both families are separated by a sur 
face which intersects the earth’s 
surface at geomagnetic latitudes not 
yet determined experimentally, but 
which according to the magnetospheric models should be about 78-80° 
in the noon meridian, and 65-7 0° at midnight, during quiet conditions. 
This limiting surface between "closed" and "open" lines is expected to 
change considerably during geomagnetic disturbances. 
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There is no convincing experimental information about the exten- 
sion of the geomagnetic tail ; at about 30 R e there is no indication in the 
data for a closing-in of the magnetopause. Theories differ widely about 
the lower and upper limits for the length of the tail. At about 20 R e , the 
tail has a diameter of roughly 40 R e . 

Taking into account the measured 
values of the magnetic field intensi- 
ty, and assuming circular shape of 
the cros s section, one obtains a total 
force of about 4 x 10 11 dynes with 
which the tail pulls on the rest of 
the magnetosphere through magnetic 
stress. This is a very considerable 
fraction of the 4 x 10^ dynes with 
which the solar wind presses against 
the entire front side of the magneto- 
sphere. In other words, a very efficient "friction" mechanism must be 
effective at the boundary of the tail. This mechanism - very likely pro- 
vided by instabilities - may also be responsible for convective plasma 
motions in the tail, necessary to explain the stretching out of the geo- 
magnetic field lines. 

The existence of a geomag- 
netic tail has opened up new view- 
points about mechanisms for the 
conveyance of low energy solar pro- 
tons from interplanetary space down 
to the earth 1 s atmosphere, particu- 
larly into the polar regions. The 
neutral sheet, on the other hand, is 
a very attractive region for the study 
of particle acceleration, and may 
well play a crucial role in the production of night side, rapid onset, 
auroras. Recent calculations have shown that particles, accelerated in 
the neutral sheet, would emerge along a very thin layer of field lines, 
with extremely small initial pitch angles, being therefore able to reach 
regions of high magnetic intensity and to precipitate into the dense layers 
of the nightside atmosphere. 

An important, but experimentally unsolved, question is that of the 
motion of the magnetospheric plasma . Or., which is equivalent in view 
of the almost infinite conductivity, the question of the motion of the 
geomagnetic field lines. These field lines are "solidly rooted" in the 
ionosphere. If the ionosphere co-rotates with the earth (which may not 
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be true at all for high geomagnetic 
latitudes), the field lines must co- 
rotate, and with them, the whole 
magnetospher ic plasma. The driving 
force for this motion, of course, must 
be provided by electric fields which 
cause char ge -independent drifts, i.e. 
bulk motions , normal to the field lines. 
The electric fields are extremely im- 
portant for the motion of moderately 
energetic particles, such as, for in- 
stance, auroral particles. 


III. MOTION OF CHARGED PARTICLES IN THE MAGNETOSPHERE 


A general study of the adiabatic motion of energetic particles in 
the earth's magnetic field is necessary for a better understanding of 
trapped particle phenomena. Let us consider a time -independent mag- 
netic field, in absence of electric fields. In this case, only two of the 
three adiabatic invariants are needed, in order to describe the longitu- 
dinal drift of a particle: the conservation of the magnetic moment, M, 
and that of the "longitudinal invariant, " J: 



ZmB 


K. E. 

= const 


( 1 ) 



const . . . ( 


2 ) 


p is the momentum, B m the field intensity at the mirror point. The in- 
tegral (2) is taken along a field line, for a complete bounce motion. In 
absence of electric fields, the particle*s energy is also conserved. In 
this case, the two invariants (1) and (2) can be replaced by two other 
invariants, which now depend on the field geometry only: 

m 1 , 1 

( B(sT 

B m = const ...(3) I = ) \| 1 - ds = const . . . ( 4 ) 

m B m 



(4) is extended along the field line 
between the two conjugate mirror 
points. For any trapping field geo- 
metry, we can assign to each point 
in space a pair of values I, B m , 
such that a particle mirroring there , 
has the value I for the integral (4), 
B m being simply the field intensity 
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at that point. As the particle drifts to other field lines, it must keep 
these values constant, i. e. , it will cover a shell of field lines which 
pass through the intersections of two given constant I and constant B m 
surfaces. 


Let us consider the geomagnetic field. Take a particle which starts 

at a given longitude, sitting on a given field line and mirroring at a value 

B m . The integral (4) computed along the field line between the two mirror 

points, has a value I. This means that when drifting through any other 

longitude, say 180° away, this particle will be bouncing along a field 

line which passes through the intersection of the corresponding I = const 

and B m = const surfaces. Now take a particle which starts on the same 

initial field line, but which mirrors at a lower value B* < B^ . Its in- 
3 mm 



INITIAL FIELD LINE 
AT 4 > 


S PLITTED SHELLS 
AT <#> + 180° 


tegral (4) will also be lower, I ? < I. After an 180° longitudinal drift, this 
second particle will be traveling along a field line which passes through 
the intersection of the surface I 1 = const and B^ = const. Only in case 
of perfect azimuthal symmetry (as in the pure dipole), will these sur- 
faces intersect exactly on the same line as that of the first particle. In 
the general case, particles starting on the same field line at a given 
longitude will populate different shells, according to their initial mirror 
point fields, or, what is equivalent, according to their initial equatorial 
pitch angles a e - arc sin(B e /B m ) (of course, all these different shells are 
tangent to each other at the initial field line). 

For the case of the real geomagnetic field (in absence of external 
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perturbations), it can be shown that 
the distance between splitted shells 
is only very small, a fraction of 1% 
of the distance of the equatorial point 
of a field line to the center of the 
earth. In other words, with a very 
good approximation, one can say that 
all particles initially on the same field 
line, will mirror on a common field 
line at any other longitude. This has 
a formidable consequence: it enables 
a two-dimensional description of the three-dimensional radiation belts, 
at least up to distances of about 5 R e . Indeed, if particles do populate 
the same shell irrespective of their initial mirror points, the omnidi- 
rectional flux of these particles will be the same on all points of the shell 
having the same B value (provided of course, no appreciable injections or 
losses occur during the drift). In order to describe omnidirectional par- 
ticle fluxes in the inner magnetosphere, we therefore need only two 
"space" parameters: the value of the magnetic field intensity at the point 
of measurement and a parameter which characterizes the (unique) shell 
which goes through that point. This is the famous L-parameter. L is a 
particular relation between I and B m which remains constant (within < 

1%) on a given field line, and, therefore, on the whole shell generated 
by particles starting on that field line. Numerically, L gives the aver- 
age distance of the equatorial points of a shell to the center of the dipole. 

But what happens in the outer magnetosphere, where the azimuthal 
symmetry is brutally removed? Particles starting on the same field line, 
say in the noon meridional plane, will now populate different shells, de- 
pending on their initial mirror points or equatorial pitch angles. For 
instance, they will cross the midnight meridian on different lines. 

Let us start with a particle mirroring at or near the equator, on a 
line in the noon meridian, close to the boundary. For this particle, I - 
0; it will drift around the earth on the equator following a constant-B 
path. This constant-B path comes considerably closer to the earth at 
the night side, because the field is weaker there (less compression), 
and we must go to lower altitudes in order to find a given B value. On 
the other hand, a particle which starts on the same field line on the noon 
meridian, but which is mirroring at high latitudes, will have a high I 
value. Under these circumstances, it can be shown that the value of I 
is not much different from the length of the field line between mirror 
points. On the midnight meridian, the particle will therefore be found 
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on a line which has nearly the same length than the initial one, i.e. 
stretching out to roughly the same equatorial distance . In summary, all 
particles initially on the same noon-line, will cross the midnight plane 
on line portions shown in the figure. Furthermore, it is easy to realize 
that particles mirroring inside that area (BB 1 ), will cross the noon me- 
ridian outside (AA ! ) the initial line. 

If this noon-line is on the boundary, 
no stably trapped particle could be 
found mirroring inside the hatched 
area in the midnight meridional 
plane. Any particle doing this would 
notbe able to complete a drift around 
the earth: it would leave the magneto- 
sphere before reaching the noon me- 
ridian. We may call this a "pseudo- 
trapped n particle (only transiently 
trapped). Notice finally that a sharp trapping boundary in the noon side 
does not result in a sharp boundary in the back side. 
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On the other hand, for a given 
field line in the midnight meridian, 
all particles mirroring anywhere on 
this line, will cross the noon merid- 
ian in an area like the one shown in 
the figure. All particles mirroring 
outside that area (BB ? ), will cross 
the midnight meridian outside (AA 1 ) 
the given line. If now there is an 
'’obstacle” behind that line (like for 
instance the neutral sheet), no sta- 
bly trapped particle could be found 
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outside the hatched area in the noon meridian. Any particle injected 
there, would be lost into the "obstacle" before reaching the midnight 
meridian; in this high latitude noon region, only pseudo-tr apped parti- 
cles can exist. 


From these considerations we deduce that the outer magnetosphere 
may be divided into the following zones: (1) a genuine trapping region; 

(2) a region in the back side, centered at the equator, which can be 
crossed by stably trapped particles, but in which they are not allowed 
to mirror; (3) a high latitude region in the front side where only pseudo- 
trapped particles can exist; (4) a region defined by the "open" field lines 
emerging from the polar caps, in which no trapping is possible at all, 
not even for one little bounce. The shape and simultaneous existence of 
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regions (Z) and (3) will depend on the position of the magnetopause and 
that of the neutral sheet, respectively. 

These considerations are of crucial importance when data on spa- 
tial flux distributions of energetic particles in the outer magnetosphere 
are being analyzed, ordered, plotted and compared. Notice finally, that 
the L -value has no physical meaning beyond L = 5 or 6; the azimuthal 
asymmetry destroys any hope for a two-dimensional description of 
omnidirectional particle fluxes in the outer magnetosphere. 

So far, we have not taken into account at all the energy of the par- 
ticles. There are, however, quite distinct regions in energy space, re- 
garding trapping conditions. For a given L- shell, there is an upper limit 
for the energy of a particle, beyond which no adiabatic motion, viz trap- 
ping, is possible. This limit is given by the energy at which the parti- 
cle's gyroradius becomes of the order of a typical scale in the field geo- 
metry. There is no clear cut between the two regions; a proof for this 
is given by the good correlation found between the n non-adiabatic M 
quantity of local vertical geomagnetic cut-off rigidity for cosmic rays, 
and the adiabatic concept of local L-value. 

On the other hand, there is a lower limit in energy, below which 
the effect of electric fields in the magnetosphere cannot be neglected. 

In this energy region, longitudinal drift velocities are comparable to 
the co- rotation velocities of the field lines. To obtain the total drift ve- 
locity, the electric drift must be added vectorially to the magnetic drift; 
anything may result. Furthermore, energy conservation is no longer 
true, and particles can be accelerated or decelerated in the electric 
field. Due to conservation of the magnetic moment (1), any increase in 
kinetic energy will be accompanied by an increase in mirror field in- 
tensity (lowering of mirror point altitude). And due to the conservation 
of J (Z), this increase will also be accompanied by a decrease in L, i. e. 

a decrease of geomagnetic latitude 
of the end points of the correspond- 
ing field line. Auroral particles be- 
long to this energy region. 

So far we have not mentioned 
at all the third adiabatic invariant, 
which expresses the constancy of the 
magnetic flux enclosed by a given 
particle shell. This invariant is im- 
portant only for time -dependent 
magnetic fields. 
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Let us say a few words about the violation of the invariants. There 
is a characteristic period of time associated to each of the three invari- 
ants. For the magnetic moment, it is the cyclotron period; for the sec- 
ond invariant, it is the bounce period, and for the flux invariant, the 
longitudinal drift period. In the geomagnetic field, these periods are in 
general orders of magnitude apart (for inner belt electrons, typically, 
~10’ 6 sec, ~10”^sec and ~10^ sec, respectively). Any time variation 
of the geomagnetic field with a time scale of the order of one of the 
characteristic periods, will lead to a violation, i.e. non-conservation, 
of the corresponding invariant. If, on the other hand, the time variation 
is much slower than a given characteristic period, it will leave the cor- 
responding invariant untouched. 

A fluctuation of the geomagnetic field with a typical time scale of, 
say, 15 minutes, will violate the conservation of the flux invariant for 
all particles having longitudinal drift periods longer than this; but it 
will not affect the conservation of M and J. If the field remains azimuth- 
ally symmetric during the perturbation, particles will always remain on 
common shells, even if they change position. The situation is quite dif- 
ferent, when there is asymmetry. In this case, everything will depend 
on where (at what longitude) a given shell-particle was surprised by the 
violation. Once the field is back to the initial value, a given shell will 
be "smeared out n in the final state. This represents a very important 
shell diffusion mechanism. Any particle which during this process is 
brought to a lower shell will increase both its energy and its mirror 
point field intensity due to the conservation of M and J. It can be shown 
that this increase will also be accompanied by an increase of the equa- 
torial pitch angle a e = arc sin (B e /B m )^ 2 (because B e increases faster 
than B m ). This means that particles get more and more confined near 
the equator, as they diffuse inwards. 

It is important to note that the ultimate physical cause for violation 
of the third invariant is given by induced electric fields, which are in- 
tense enough to make particles drift out of their "home shell". Energy 
changes are precisely betatron-type accelerations in these fields. 

Let us finally consider an extremely short perturbation, like the 
elastic scattering of a trapped particle. In this case, the first and the 
second invariants will be violated and the particle will change its mir- 
ror point. In the case of symmetric field, the third invariant remains 
almost unchanged: the particle will remain on the same shell within the 
accuracy of one cyclotron radius. In an asymmetric field, however, we 
already have seen that a change in the mirror point of a particle leads 
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to a change of the particle's shell. Particles diffusing in pitch angle 
will therefore diffuse across shells too, in an asymmetric field. This 
type of radial drift, however, will not be accompanied by any change in 
energy (this time, there is no electric field around to accomplish the 
acceleration) . 

IV. GEOMAGNETICALLY TRAPPED RADIATION 

A great deal of information about the shape and the structure of 
the magnetosphere comes from the study of energetic particle distribu- 
tions in space and time. We will discuss here only some very general 




features regarding radiation belts. 
First of all, if one estimates the 
ratio of particle kinetic energy den- 
sity to magnetic field energy density 
(a quantity usually called /?, indic- 
ative of the score in the game "par- 
ticles vs. fields"), one obtains very 
high values, reaching 0. 1 beyond 
L ^ 3 (labor atory plasma experiments 
never achieved more than 10“ 3 ). 

This indicates that the earth's radia- 
tion belts are probably loaded up 
with the maximum flux of particles 
compatible with plasma stability. 

We first turn to the spatial dis- 
tribution of energetic particle fluxes 
in the magnetosphere. It is quite well 
established that below L ~ 6, omni- 
directional flux contours give a con- 
sistent picture in the two dimensional 
B - L space, indicating that shell 
splitting due to azimuthal asymmetry 
is negligible. An over -all qualitative 
picture of omnidirectional particle 
fluxes, near the equatorial plane, is 
shown in the sketch. The "inner" ra- 
diation belt with peak fluxes at about 
L = 1.6, is very stable, both during 
geomagnetic storms and throughout 
the solar cycle. The outer belt shows 
considerable time variations, mainly 
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correlated with geomagnetic activity. Notice the characteristic "slot" 
in the electron f luxe s around L = 2. 5 - 3, more pronounced for higher 
energies. The flux of high energy protons decays very rapidly with L, 
likely due to the sharp decrease with increasing radial distance of the 
upper limit in energy for adiabatic trapping of protons. 

Energy spectra of protons are now quite well known. In the inner 
belt, measurements extend to almost 1 Gev. In the outer zone, where 
high energy protons are absent, exponential spectra are found, with 
e -folding energies ranging from tens of Kev at high L values, to 400 
Kev at L = 3 . A very clear L - dependence of these spectra is thus es- 
tablished, indicating that protons drift in L, hardening their spectra 
towards lower L shells, as dictated by the conservation of M and J. In- 
tense fluxes of 5 Kev protons have recently been found at low L values. 
They may play an important role for the ring current during geomag- 
netic storms. Energy spectra of natural electrons in the inner belt are 
not so well known; contamination by artificially injected electrons is 
still interfering the measurements. In the outer zone, spectra may ten- 
tatively be fitted by an exponential form with e -folding energies ranging 
from about 170 Kev at high L values to about 340 Kev at L = 4. Again, 
there is a trend for the spectrum to harden towards lower L values. 

Beyond L = 6, the distortion of the magnetosphere calls for a three- 
dimensional description of the radiation belt. In the equatorial plane, 
isointensity contours consistently follow constant B rings, coming closer 
to the earth at the night side. Most of the information in this outer region 
comes from >40 Kev electron measurements. Typical fluxes here are 
of the order of 10 7 cm " 2 sec " 1 . 

The "core 11 of trapped particles shows in general a well defined 
boundary. In the equatorial plane, 
and near the noon meridian, this 
trapping boundary is very close to, 
or coincident with the magnetopause; 
it shows up as a sharp drop in in- 10 6 _ 
tensity by several orders of magni- 
tude, over a rather short distance. 

As one goes towards the dawn or 
evening sides, the trapping boundary 
and the magnetopause get disconnec- 
ted, leaving between them a region 
where strongly agitated fluxes are 
observed; typical intensities of >40 
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Kev electrons are here of the order 

4 - 2 - 1 

of 10 cm sec . These electrons 
have much softer spectra than in the 
trapping region. They may be pseudo- 
trapped particles, injected or accel- 
erated in the day side, bouncing 
along closed fields lines, but unable 
to complete a drift around the earth 
(page 11). In the night side, always 
on the equatorial plane, there is a 
much smoother and less marked 
transition from the high fluxes of 
stably trapped electrons into the un T 
stable region, occuring in general at about 8 R e (see discussion page 11). 

This region of strongly varying fluxes 
of soft electrons extends out to about 
15 R e into the tail, in the sun -earth 
direction. 




EQUATORIAL 

PLANE u Tr ' * nr ' rri i\\\\\\m The w hole picture changes 

greatly at higher geomagnetic lati- 
tudes. The trapping boundary moves 
inwards, as expected in view of the 
magnetic field configuration. On 

the day side, the higher energy {>280 Kev) electron fluxes seem to ter- 
minate at lower latitudes than the 
> 40 Kev flux. Going to the night side, 
the region of strongly varying soft 
electron fluxes stretches out like a 
"skirt" along the sides of the mag- 
netosphere. Near the midnight merid- 
ian, at these higher latitudes, this re- 
gion disappears completely: the stably 
trapped electron flux drops abruptly 
to the background. This clearly in- 
dicates that the "cusp" of soft elec- 
tron fluxes stretching out into the tail, 
mentioned in the preceding paragraph, is strictly confined to low lati- 
tudes. It is probably related to the presence of the neutral sheet. 

A very important fact is the appearance of electron "islands" in 
the back side of the magnetosphere. These are sharp increases of 
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electron fluxes, up to 1 0 7 cm " 2 sec " 1 
for > 40 Kev, followed by an expo- 
nential decay. They are clearly time- 
dependent phenomena. Similar events 
were observed near the magneto - 
pause and in the transition region. 

Their frequency of appearance is 
positively correlated with K p0 

Recent measurements of very 
low energy electrons (0. 1-10 Kev) 
show a distribution profile of high 
fluxe s (10® - 10 9 cm ” 2 sec -1 Kev” ^ ) , 
strongly dependent on geomagnetic 
activity, which in general seems to extend from low latitudes in the back 
side towards high latitudes at the dawn and evening sides. Turning finally 
to extremely low energy electrons, in the range of 100 ev, no clear 
boundary is seen at all, until the shock is reached. This seems to in- 
dicate that the magnetopause is quite transparent to these particles. 

The analysis of time variations of the radiation belts is a highly 
complicated business. We can only refer here to some of the salient fea- 
tures. Only the core of trapped particles will be considered. First of 
all, we may mention the following types of time variations : the 11 year 
variation; transient variations associated with storms or high-K p ; de- 
cays after artificial injections, and short term periodic variations. The 
study and the interpretation of these variations is complicated by the 
fact that it is very difficult to observe the flux and energy spectrum at 
a fixed point in the magnetosphere throughout a typical time scale of the 
variation; hence it often is very difficult to determine whether one is 
dealing with a variation in flux or in energy spectrum, whether one has 
just a shift of the whole spatial distribution of the belt, or whether sev- 
eral of these occur simultaneously. 

Regarding the 11 year variation, it seems an established fact that 
the inner belt is pretty constant throughout the solar cycle. The outer 
belt particles show so many and intense short term variations, that it 
is difficult to establish a "base line" which then could be followed 
throughout the solar cycle. In general, a gradual decay of intensity 
during declining solar activity was detected. In any case, it seems clear 
that what keeps the outer belt alive, is the solar wind, or more pre- 
cisely, the solar wind fluctuations. Turning them off, very likely 
the outer radiation belt would disappear. 
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As to the short term variations associated with increases in the 
geomagnetic activity, one may say in summary that they are more in- 
tense at higher L. values, that they are felt more by electrons than by 
protons, that variations of high energy particles are delayed in time up 

to several days with respect to the 
lower energy ones, and that certain 
common patterns in time seem to be 
followed by electrons and protons, 
respectively (see sketch). One common 
feature is a sudden change, followed 
by a gradual recovery. The relax- 
ation times of these recoveries pro- 
vide important information on par- 
ticle lifetimes. Typical maximum 
increases for >40 Kev electron 
fluxes are given by factors of the 
order of 100. Another important in- 
formation came from studies of in- 
ward motion of intensity profiles of 
outer belt electrons, after a big 
storm. Assuming that this repre- 
sents a diffusion in L, one can de- 
rive values for the diffusion M veloc- 
ity n dL/dt. 


Artificial injections of elec- 
trons by means of nuclear explosions 
have provided important information 
about particle lifetimes. The com- 
monly observed feature is a rather 
fast decay of the initially strongly 
antisotropic distribution of particles 
into the "normal mode" correspond- 
ing to each L shell; this normal mode 
decays roughly exponentially, with a 
characteristic lifetime which not only 
depends on L, but also depends on the 
mirror point field intensity ( or equa- 
torial pitch angle) of a given parti- 
cle. Typical lifetimes for 2 Mev 
electrons are sketched in the figure, 
as a function of L. 
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A very important type of periodic time variations of energetic 
electron fluxes was found recently, occuring after abig storm, in asso- 
ciation with synchronous variations of the magnetic field. Harmonic 
analysis of the flux variations revealed a strong predominance of a fre- 
quency, which roughly coincided with the mean longitudinal drift frequency 
of this group of electrons. This suggests the action of a resonant-type 
acceleration in a periodically varying field. Such a process, of course, 
is possible only if the induced electric field is azimuthally asymmetric, 
in order to give chance to a bunch of particles to consistently see a 
stronger field during the accelerating phase, every time they drift 
around the earth. This asymmetry may be provided by the natural dis- 
tortion of the magnetosphere. 

In regions of steep flux gradients, as at the lower edge of the inner 
radiation belt, this process may lead to considerable effects, even if 
there is only a very small energy increase of the resonant particles: as 

they get accelerated, they will drift inwards in L (conservation of M and 
J), i.e. to substantially lower flux regions, where they may stand out 
considerably over the background, as a monoener getic group. This effect 
was observed recently in the form of sudden appearance and gradual de- 
cay of a 1.3 Mev electron peak, at an L shell of about 1.15. Equatorial 
magnetograms during the first hours showed marked recurrent variations, 
with a period in resonance with the 1.3 Mev electron drift around the earth. 

The last part of our discussion will deal with radiation belt dy- 
namic s , i.e. with the sequence of processes injection -►storage — lo"ss. 

Let us consider a given class of particles. Under storage, we mean the 
history of the particle from the instant it becomes a trapped particle 
(injection), to the moment it disappears from the scene (loss). During 
this time, the particle is trapped in the geomagnetic field, subject to 
pitch angle scattering, L- shell diffusion and acceleration. Physical 
mechanisms governing these processes entirely determine all properties 
of the radiation belts. 

Let us first start with the loss mechanisms . The ultimate sink, at 
least for lower L shells, is always the atmosphere. A particle, trapped 
on a given shell, may leave it without precipitating into the atmosphere. 

In that case, it just was transfered to another neighboring shell. But 
some day, somewhere, it will die in the upper atmosphere. In the inner 
belt, L-transfer processes seem to be very weak: it is more likely for 
a particle to get the n right n kick in pitch angle through a scattering 
process, ending up in the dense atmosphere. In the outer belt, during 
quiet conditions, L-diffusion (towards lower L shells) seems to be the 


19 


dominating process; losses into the atmosphere are enhanced only during 
geomagnetic disturbances, which very likely make the pitch angle scat- 
tering mechanism temporarily overwhelm the L-transfer process. For 
low L- shells, most of the precipitation into the atmosphere occurs in the 
M South American" or "South Atlantic" Anomaly, where a given particle 
shell has its closest approach to the earth's surface (mainly because of 
the eccentricity of the geomagnetic dipole). Electrons, drifting from west 
to east, attain lowest mirror point altitudes in this Anomaly; those getting 
below, say, 100 km will be wiped out from the radiation belt by energy 
loss in the dense atmosphere. This leads to a region in B-L space with 
depleted electron fluxes, east of the Anomaly. However, this region is 
again replenished to some extend with electrons which diffuse into it by 
pitch angle scattering. These electrons will precipitate into the atmosphere 
the next time they drift over South America or the South Atlantic. A 
permanent flux of X-rays from electron Br ems str ahlung has indeed been 
detected at balloon altitudes in that area. 

There are at least two types of pitch angle (or mirror point) dif- 
fusion mechanisms. One is elastic, multiple Coulomb scattering with 
air atoms. This process is only effective for electrons (for protons, 
slowing down by ionization loss and, eventually, charge exchange, 
are dominant). It describes very well all experimental observations for 
trapped electrons at very low altitudes (L^l.Z). For higher L values, 
its effects must become negligible, because of the extremely low atmos- 
pheric density. However, there is definite experimental evidence that a 
second efficient pitch angle scattering mechanism exists, and even in- 
creases with higher E values. The existence of the "slot" in the electron 
distribution at L - 3 is one indication in favor. The short lifetimes of 
electrons beyond L = 3 is another evidence. More support comes from 
experimental results about the efficient replenishment of electrons east 
of the South American Anomaly, in all B-L regions which plunge below 
sea level in the South Atlantic, i.e. which get depleted there. Calcu- 
lations have shown that Coulomb scattering must refill these regions 
partially, a few degrees east of the center of the Anomaly; however, 
satellite measurements indicate that refilling c ontinue s at longitude s where 
Coulomb scattering is completely inefficient. A possible pitch angle 
scattering mechanism is given by resonant interactions of electrons 
with electromagnetic waves, such as whistlers. This is particularly at- 
tractive, for it predicts a maximum efficiency near L = 3, where the 
slot occurs. However, theoretical predictions so far show no agreement 
with experimental data on pitch angle distributions. Another possibility 
would be interaction with radiowave noise, or with hydromagnetic waves. 
Protons would not be affected by these mechanisms. 
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The other extremely important mechanism is that of L -diffusion. 

We already have mentioned, that the remarkable L-dependence of the 
shape of proton spectra is a strong indication that we are dealing with 
one and the same population of particles, which move inwards in radial 
distance. A mechanism which explains quantitatively this diffusion is 
provided by the violation of the third invariant (page 13) during periods 
of enhanced geomagnetic activity. This, together with shell splitting and 
conservation of M and J, leads to a gradual diffusion of protons towards 
lower shells, increasing their energy, and increasing their equatorial 
pitch angles. It is a rather slow process; it takes years for a proton to 
diffuse to low L values - this may explain the stability of the inner 
belt with respect to the 11 year cycle. 

Electrons very likely undergo similar shell diffusion processes. In 
this case, however, the efficient action of pitch angle scattering mech- 
anisms blurrs the picture considerably. As explained on page 14, any 
diffusion in pitch angle in a strongly asymmetric field will be accompa- 
nied by shell diffusion, even in the absence of geomagnetic perturbations. 
In this process, however, energy will be conserved. There is no theory 
yet, for the electron belts. 

In addition to the betatron-type acceleration of a particle which 
moves across E- shells during third invariant violations, there must be 
other acceleration mechanisms. One already mentioned is the resonance - 
type mechanism acting during periodic variations of the earth 1 s field. On 
the other hand, the characteristic time variations of the outer belt during 
Kp increases, may point to a third type of a locally acting acceleration 
mechanism. 

As to injection mechanisms, most of the trapped particles ulti- 
mately may come from the solar wind. They may be brought into the 
magnetosphere at very low energies through neutral points, or instabilities 
at the magnetopause, being then energized as they drift towards lower 
shells. On the other hand, part of the belt population may be magneto - 
spheric plasma, locally accelerated. Finally, cosmic ray albedo neutron 
decay certainly is a necessary, but not a sufficient, source. It may be 
sufficient for very high energy protons. It fails, however, to explain the 
observed fluxes and spectra, and solar cycle variations of intermediate 
and low energy particles. There is really only one well known source of 
geomagnetically trapped particles: the bomb. 
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